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In |LY] a differential Harnack inequality was proved for solutions to the heat 
equation on a Riemannian manifold. Inspired by this result, Hamilton first 
proved trace and matrix Harnack inequalities for the Ricci flow on compact 



surfaces [H0| and then vastly generalized his own result to all higher dimensions 



for complete solutions of the Ricci flow with nonnegative curvature operator 



[ H2]. Soon afterwards, a matrix Harnack inequality for the Kahler- Ricci flow 
under the assumption of nonnegative bisectional curvature was proved by Huai- 
DongCao (C). In this paper, following a suggestion of Richard Hamilton, we give 
an alternate proof of the matrix Harnack inequality for the Ricci flow originally 
proved by himj^ In | H2|, Hamilton proved that if {M",g{t)) is a complete 



solution to the Ricci flow with nonnegative curvature operator, then 
dR R 



dt t 



2V^RV' + 2i?,jW^ > (1) 



for any vector field V. In particular, if in addition the Ricci curvature is positive, 
then 

This trace Harnack inequality is a corollary of Hamilton's matrix Harnack in- 
equality, which we now recall (we follow Hamilton's sign convention for the 
curvature tensor so that Riju — gkhR'iji)- Define tensors P and M by 



Mab 4= ARab — -DaDbR + 2RacbdRcd ~ RacRbc + '^^ab 



•Research partially supported bvJNSF grant DMS-9971891 



^Yet another proof is given in [CC|. Hamilton's motivation for considering the alternate 
proof given in the present paper is that it may possibly falicitate extending his Harnack 
ineg ualitv to the variable signed curvature case, especially in dimension three. See the end of 
JHSl for some reasons for studying the problem of extending the Harnack inequality. 
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and consider the Harnack quadratic 

Z ({/, W) # M,^,WaW^ + 2P,b,C/,6W, + RabcdUabUcd, 

where U is any 2-form and W is any l-form. Hamihon's matrix Harnack es- 
timate says that if {M^^,g{t)) is a complete solution to the Ricci flow with 
nonnegative curvature operator, then 

Z [U, W)>Q 

for all U and W. Taking U — V AW and summing over W in an orthonormal 
frame, one obtains (|]). 

The idea of giving an alternate proof of it is as follows. Given any vector 
W, consider the 2-form U minimizing the quadratic form. This leads to the 
consideration of a symmetric 2-tensor Zab which has a nice evolution equation. 
In particular, Zab is essentially a supersolution to the heat equation so that 
the maximum principle implies that it is positive definite, which is equivalent 
to the matrix Harnack inequality. Tracing this 2-tensor yields a trace Harnack 
inequality: Z = Zaa > 0. 

Now we proceed with the details of the proof. Assume that the curvature 
operator is positive and let Sabcd denote the inverse of Rabcd so that 

^ ^ Sabef Ref cd — ^abcdj 

where / : A^M h^M is the identity. First we rewrite the matrix Harnack 
quadratic by completing the square: 

Z {U, W) = MabWaWb - S.jklPjpWpPklqWg 

+ Rabcd {Uab + SabijPijpWp) (Ucd + ScdklPklqWq) . 

Since Rm > 0, the above expression is minimized when 

Uab + SabijPijpWp = 0, 

in which case 



Let 



Z {U, W) = {Mab - S^jklPjaPklb) WaWb. 



Zab — ^'lab ^ SijklPijaPklb, 



which is a symmetric 2-tensor. 

The main computation is that in an orthonormal moving frame, the evolution 
of Zab is given as follows. Following the idea in | H2 1 , as long as it is nonnegative- 
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definite, we may write the Harnack quadratic as a sum of squares: 

^abcd — J^ab^cd 
N 

^abc — 2^ ^ab^c 
N 

Mab = Y.^-^b 



^b ■ 



N 



Here, given a point x G M, {Y^ is some orthogonal set of vectors 

in A^M^;© A'^M^ and m = dim A^M^, + dim A^M^ = nll^. 

Main Theorem. // Rm > 0, then the Harnack quantity Zab satisfies the 
evolution equation: 

^ ~ ^ j Zab — '^SmntuKavtuKbvmn + L^^^L^^^ — —Zab, (2) 



wher 



^avtu — Sijrs-Rija^ vRrstu ^v-Rtua RtuawRvw ~t~ ^^Rtuav 
— ^^id ^jd 'Jijkl^kla + J^ah^h ~ ^ah^h ■ 



Note that applying the maximum principle for symmetric 2-tensors in section 



7 of I HI I yields Zab > 0, which is equivalent to Hamilton's matrix Harnack 
inequality. Here we used the fact that if Rm > 0, then as long as Zab is 
nonnegative definite, i^*^ is well-defined. 

The bulk of the rest of this paper is devoted to proving the theorem. First re- 
call Hamilton's computations for the evolution equations of the Harnack quadratic 
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Lemma 1 In an evolving local orthonormal frame field, we have 
1. 

{Df — A)Rabcd — 2{Babcd ^ Babdc + Bacbd ~ Badbc) , (3) 

where Babcd = 

2. 

[Dt — A)Pabc = —2RdeDdRabce + 2RadbePdec + ^RadcePdbe + "^RbdcePade ■ 

(4) 



3. 



{Dt - A)Mab = 2R,d[DcPdab + D,Pdba] + 2RacbdM,d (5) 
+ 2Pa.dPbcd - APacdPbdc + 2R,A,Radbe ^ ^ Rab ■ 
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We use this lemma in the computations below. We start by computing: 

— Sijkl - A^ PijaPklb - SijklPija " Pfejfe 
+ 'iSijkl^ pPija^ pPklb + '^^pSijkl'^pPijaPklb + '^'^ pSijklPija^ pPklb- 

Thus, besides the evolution equations for Mab and Pabo we need the evolution 
equation for Sijki , which is given by: 

^^ijvs^mntu^klpq ~t~ Sijran^klrs^pqtu) ^ yPrstu^ vPmnpq 
— '^^ijmn^klpqi^-^mnpq ^mnqp ~\~ -^mpnq -^mqnp^ 

{^ijrs^mntu^klpq ~t~ ^ijmn^klrs^pqtu) ^vPrstu^vPmnpq- 

Substituting in those evolution equations, we obtain: 

d \ 

— — A j Zab = 2Rcd[DcPdab + DcPdba] + '^Racbd.Mcd 

+ '^PacdPbcd — ^PacdPbdc + ^PcdRcePadbe — Rab 

{Sijrs^mntuSklpq SijmnSklrsSpqtu) ^ v Rr stu^ v Rmnpq 
~\~'2SijrfYi'[iSklpq{^Bjyijipq -^mnqp ^ Pmpnq Prnqnp) 

— "^SijklPklb { — RdeDdRijae + R-idjePdea + RidaePdje + RjdaePide) 

— "^SijklPija { — RdeDdRklbe + RkdlePdeb + RkdbePdle + RldbePkde) 

~\~ ^Sijkl^ pPija^ pPklb '^Sijmn^klpq^rRmnpq rPijaPklb ~t~ Pija^ rPklb^ 

where Babcd = RaebfRcedf ■ Taking into account symmetries and combining 
some like terms, we have 

d 

— — A I Zab = "^RacbdZc, 



PijaPklb 



dt 



^SijjYifiS]^lpqS/ j-RjYinpq r PijaPklb ~\~ Pija^ rPklh) 
~\~ {^ijrs^mntu^klpq ~l~ ^ijmn^klrs^pqtu) PijaPklb^ v-^rstv^ v-^^^npq 
+ "^RacbdSijklPijcPkld + "^PacdPbcd — ^PacdPbdc 

-\- '2Sij'jnnSklpq{Pmnpq Pmnqp H" Prnpnq Pmqnp^ PijaPklb 
~ 2 {SdeklPdiej H" SijdeRdkel) PijaPklb 

'^^ijklP'idaePklbPdje ~ ^Sij^^^lRj^dbePijaPdle 

+ "^RcdlDcPdab + DcPdba] + ^PcdRceRadbe — P-ab 
+ "^SijklPklbPdeP^dP-ijae + Sijkl Pi jaPdeDdP-klbe- 
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Separating the terms quadratic in Pabc from the rest and noticing that most of 
the rest is a square, we obtain: 

|-^)^- 

= '^Racbd^cd 

~l~ ^^mntu {^ijrs-Pija^v-^rstu ^v^tua ~l~ -^uaw^vw^ 

+ '^RacbdS ijkl PijcPkld + '^PacdPhcd — ^PacdPhdc — "^SijklRidaePklbPdje 

^^ijklPkdbePijaPdle ~t~ '^SijmriSklpqi^^mnpq ^mnqp ~l~ -Pmpnq ^mqnp^PijaPklb 

~ 2 {S(iQfi;iR(iiQj + SijdeRdfi^Qi) PijaPklb ~ Rab- 

We simplify the terms that are quadratic in Pabc • 

A =2RacbdSijklPijcPkld + '^PacdPbcd — ^PacdPbdc — "^SijklRidaePklbPdje 
~h '^SijmnSklpqi^Pmnpq Pmnqp H" P nipnq Pmqnp) P'ijaPklb 

— 4:SijklRkdbePijaPdle — 2 (SdeklRdiej + SijdeRdkel) PijaPklb 

— '^Racbd^ijklPijcPkld H" '^PacdPbcd ^PacdPbdc 
^^ijklRidaePklbPdje ^^ijklRkdbePijaPdle H" PcdaPcdb 

~\~ 2Sijmn^klpq{,Pmpnq Pmqnp^ PijaPklb PklaPklb PijaPijb 
= —'^PacdPbdc + '^RacbdSijklPijcPkld — ijklRidaePklbPdje 

^^ijklRidbePklaPdje H" '^RiepjRjeqfSijj^lPfQlaSpquitiPmnb 

^Rieqf Rjepf Sijkl Pkla S pqnin Prnnb 

— '^PacdPbdc ~\~ '^Racbd^ijkl PijcPkld ^SijkiRidaePklbPdje 
^^ijklRidbePklaPdje H" ^RiepfRjeqf ^ijklPkla^pqmnPmnb- 

In the above computation we used the foUowing identities. First, the quadratic 
B in curvature can be rewritten using 

Rmnpq ~ RmnefRpqef ~ {Rmenf Rnemf^ (.Rpegf Rqepf^ 
~ Pmnpq P-mnqp Pnmpq H" Pnmqp 
~ 2 {Pmnpq Pmnqp^ • 

and 

'^SijfnnSklpqi^Pmpnq Pmqnp) PijaPklb 

— '^SijmnSklpqi^Rmepf Rneqf Rmeq f Rnepf) PijaPklb 

— '^SijklS'tnnpqRiepf Rjeqf PklaPmnb '^Sijj^iSfjinpqRieqf Rjepf PklaPmnb- 

Second, we have the identities 

'^SdeklRdiej Pija ~ Pkla 
^SijdeRdkelPklb = Pijbi 
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where the first follows from 



^dekl {-^diej ^djei) -^ija ^dekl-^deij-^ija ^kla 

which is obtained from 

and the second is proved similarly. Third, we have the identity 

PcdaPcdb = {Padc + Pcad) (Adc + Pcbd) 
= 2 {PacdPbcd " PacdPbdc) ■ 

Now define Eija = SijkiPkia to simplify notation further. Then 

— "^PacdPbdc ~^ '^PacbdSijklPijcPkld ^Pddae^ijbPdje 
^Pidbe^ijaPdje ~t~ ^PiepfPjeqf^ija^pqb • 

Recall that 

^abcd — J^ab^cd 
N 

Pabc = ^ yJ^X^ 

N 

JV 

Thus 

A = —'^PacdPbdc + '^RacbdSijklP'ijcPkld + ^YpfYalX(^ ^pqb 
^^id '■be ^ jd ^ija + 'i^ id ^ pj ^ jd ' qf ^ija^pqb 

— 4 ['^id J^jd ^ija + 'ah^h ) Vpf'^qf ^pqb + ^6e ) ~ ^'ah^h 'be ^e 
^PacdPbdc ~\~ ^Pacbd^ijklPijcPkld 

= 4 {Y^Y^^E^j^ + yj^Xf ) {Y^'jY^jE,,, + n^Xf ) - 2R^^,dZed 

- 2RacbdMcd - 2PacdPbdc- 

Note that the terms Y^Yj^Eija and Y^Y^jEpgb are both antisymmetric in N 
and M. 

Now the terms on the last line may be rewritten as 

— 2RacbdMcd — 2PacdPbdc 

^x^Nt^N yM Y^vM Y^ T3 P 

— ~^'^ac'bd^c ^d ~^'^ac^d '^bd^c ^^acd^bdc 

— (vN yM I yN\ t\rN yM I yN\ 

— ~\'^ac^c + 'ac ) {'bd^d + 'bd ^d ) 
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so that 

A = 4 {Y^YJ^E,,^ + Y^iXff) {y;}Y^JE,,, + Y,^X^) - 2R^^,,Z^ 

- {YZxr + F„f ) (yjxf + Y,^X^) 

— K'^^id ^jd ^ija + J^ah^h ~ ^ah -^h ) [^^pf^gf ^Pqb + ^be -^e ~ ^be ) 
~ '^RacbdZcdt 

where we used 

V-W (vN yM , T^M yN\ _ 1 (t^N yM I vM yN\ fvN yM , t^M yN\ 

^ah^h ['^be + '^be ) - 2 K^ah^h + ^ah^h ) [^be + ^be j 



and 



-'id ^jd ^ija [^be + ^be ) — ^ 



(since Y^'^Xi^ + F^f is symmetric in N and M). 
Therefore the evolution equation for Zab is: 

— '^Smntu i^SijrsPija^ vRrstu ~ ^v^tua ~(~ RtuawRyw) 
^ klpq Pklb^ V ^pqmn vPrnnh H" RmnhxRvx^ 

where Eija = SijkiPua- We may rewrite this as 

— '^S'mntu ^^ijrs^ija^ v^rstu ^ v^tua H" -^tuaw-^vw H" 2^"^* 

^ i^^klpqPklb^ '0 R'pqrnn ^v-Prnnb H" ^mubxRvx H" '^Rrn/nbx 

~ ^^ijrs^ija^ v^rsbv ^ v^hva H" ^bvaw^vw ~l~ '^■^bvav 
~ I ^klpqPklh^ v-Rpqav ^ v^avh ~^ ^avbx-^vx H~ "^^avby 
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Now 



^ vPbva ~t~ Rbvaw Rvw i^^^hvav 



-V„ {VbR'va — ^vRba) + RbvawRvw + 'j^^b 
^Rba ^ b^ V Rva RbwRwa ~t~ RbvawRvw 
~t~ RbvawRvw ^" 2t^^^ 

= Ai?f,a — -^b^aR + "^RbvaivRviv ~~ RacRbc + 



and 



S/yRr 



ibv 



^ rRvt 



bv 



VsRr 



jbv 



-VrRsb + Vsi?r 



if)- 



Hence 



9t 



2Sn 



A Z 



a6 



m7itu I ^ijrs^ij 



^vPi 



tua ~t~ RtuawRvw 



2t 



Rt-u 



^ I ^klpqRklb^ vRpqmn ^v-Rninb ^ RmnbxRvx ~t~ ^^Rmnbv 



+ {"^^id^jd ^ija + Y^X^^ ~ Y^i^Xj^) {2Y^fY^ Epqb + F^^X 

^ ~ i^SijrsPijaPrsb + Mab) ~ — S klpqPklbPpqa + Mab) 



That is 
di 



^vPu 



X SklpqPklh^ vRp 



V P 



2t 



:R-n 



ibv 



2t 



:Rt 



This completes the proof of the main theorem. 

Now it is interesting to observe the form of the resuhing evolution equation 
when we trace the above computation. Define 

Z=Zaa = l (^AR + 2 \Rcf + - S,,klP^JaPkla. 

This trace Harnack quantity is slightly different than (|l|) ; however its positivity 
is of course still a corollary of Hamilton's matrix Harnack estimate. Tracing 



equation (|^), we have 



(^-Q^ — Z — 2SijklKmnijKmnkl + (-^^^^) ~ -Z > ^^Z- 

Note however that an appHcation of the maximum principle does not yield 
a direct proof of the trace Harnack inequality. This is because in order for 
, which appears on the rhs, to be well-defined, we need the matrix Harnack 
quantity to be nonnegative! 

Acknowledgment. We would like to express our gratitude to Richard 
Hamilton for suggesting to us the method of proof given in this paper. 
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